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^ ^ We present a microscopic theory of ground-state spectral function of bilayer quantum Hall systems 

■^Uj . that includes interactions between Hartree-Fock quasiparticles and quantum fluctuations of the 

' order parameter fleld. The collective modes in these systems are properly described only when 

fluctuations in direct and exchange particle-hole channels are taken into account. Using an auxiliary 
field functional integral approach, we present a generalization of the random phase approximation 
for quasiparticle self-energy which captures fluctuations in both channels. We discuss its relationship 
to diagrammatic perturbation theory and an adiabatic approximation. We present simple analytical 
results for the quasiparticle self-energy and the renormalized order parameter that follow from this 



C/3 ' theory. 
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I. INTRODUCTION 



The physics of a two-dimensional (2D) electron gas in a perpendicular magnetic field is unique in many respects. 
Bilayer electron systems consist of a pair of 2D electron gases separated by a distance d{d ~ 100 A) which is comparable 
^ , to the typical distance between electrons within one layer. In the presence of a strong magnetic field perpendicular 
Q ■ to the layers, the kinetic energy of electrons is quenched and the physical properties of the system are determined by 
O , electron-electron interactions alone, jBilayer quantum Hall systems at Landau level filling factor j-u = 1 exhibit a rich 
variety of broken symmetry statesETQ and have been extensively studied over the past decadeErii3. For filling factor 
\ V = \, the ground state is fully spin polarized and the spin excitations are gapped because of Zeeman coupling of 
^ . electron spins to the external magnetic field. The only remaining dynamical degrees of freedom are the discrete layer 
' index and the intra-Landau level orbitals. We model the layer index using a pseudospin label where "up" denotes 
. the symmetric bilayer state (S) and "down" denotes the antisymmetric bilayer state (AS). Note that this choice of 
quantization axis is different from the usual convention where "up" denotes a state localized in the top layer. At total 
filling factor v = 1, because of strong interlayer correlations, the ground state of the system exhibits spauinneous 
interlayer phase coherence (easy-plane ferromagnetism in the pseudospin language) and is incompressiblcE'tll. At a 
^ \ critical layer separation d^r the system undergoes a phase transitijOH-jfrom this incompressible quantum Hall state 
^ . with pseudospin ferromagnetism to a disordered compressible statelD'tll, possibly with other more exotic intervening 
' statesB. 

^ , In this paper we present a theory which describes effect of quantum fiuctuatims on the Hartree-Fock quasiparticles 
I of the V = \ bilayer. Some of the results presented here were published earlierEj. Our theory is based on an approx- 
' imate expression for the quasiparticle self-energy which includes the effect of quantum fluctuations and generalizes 

the random phase approximation (RPA) to cases when both direct and exchange fluctuations are important. This 
self-energy modifies the spectral function of the mean-field quasiparticles, which in turn renormalizes the physical 
properties of the system. The most important macroscopic parameter characterizing the phase-coherent state is the 
pseudospin polarization. We use pseudospin polarization normalized to its mean-field value as a dimcnsionless order 
parameter m^. This order parameter is related to quasiparticle spectral weights integrated up to the Fermi energy, 
$^ ' ruz = {ns — nAs), where denote occupation numbers for the symmetric and antisymmetric states. Therefore 
5^ the order parameter is sensitive to changes in the quasiparticle spectral function. In the mean-field approximation 
the quasiparticles have sharply defined energies and the spectral functions are (5-functions at these energies. In the 
mean-field ground state, all electrons occupy the symmetric state in the lowest Landau level, which corresponds to 
all pseudospins pointing along the positive z-axis, = 1. We will find that when we include the effect of collective 
fluctuations, the spectral functions develop a continuum piece and the spectral weight at the (5-function peak is re- 
duced. We will calculate the renormalized order parameter by evaluating the fluctuation correction to the mean-field 
spectral function. 

The generalized random phase approximation for quasiparticle self-energy presented in this paper can be usefully 
adopted to study the effect of quantum fiuctuations on mean-field quasiparticles whenever the collective- mode behavior 
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is determined by fluctuations in more than one particle-hole channels; for example, to analyze the effect of fluctuations 
on Wigner crystal mean-field quasiparticles at low filling factors. Traditionally random phase approximations have 
been used either to take into account the screening of Coulomb interactiorJla where onlyj-duject-channel interactions 
determine the physics, or the effect of spin-waves on a ferromagnetic mean-field statalj'LJ where the low-energy 
physics is solely determined by exchange interactions. The generalization which we present here systematically takes 
into account fluctuations in both channels. We study the bilayer quantum Hall system as a test case because its 
translation invariance combined with the absence of dispersion in the Landau level bands permits many elements 
of this calculation to be performed analytically, and because the effect of quantum fluctuations on the mean-field 
ground state can be tuned over a large range simply by changing the layer separation. In the d — limit, since 
the interaction Hamiltonian is pseudospin invariant, the mean-field ground state is exact. At finite d, because of 
the difference between the intralayer Coulomb interaction Va{<i} = 27re^/g and the interlayer Coulomb interaction 
Ve{(1) — VA{q)e~'^'^, the symmetry of the interaction Hamiltonian is reduced from SU(2) to U(l), and the order 
parameter ruz does not commute with the interaction. Therefore the pseudospin-polarized mean-field state is not an 
eigenstate of the Hamiltonian and in particular, it is not the ground state. The exact ground state of this system 
incorporates quantum fluctuations in the quantum Hall ferromagnet's order parameter field. 

The paper is organized as follows. In section ^ we describe the functional integral approach to this problem. It is 
known that the collective modes io, these systems are properly described only when both electrostatic and exchange 
fluctuations aralaken into accountSfl. No standard Hubbard Stratonovich (HS) transformation can treat both effects 
simultaneouslyEj. Section |l| also summarizes the generalized Hubbard- Stratonovich which we employ and establishes 
some of the notation which we will use. In this section we present the generalized random phase approximation (GRPA) 
for the quasiparticle self-energy, which takes into account fluctuations in both channels and is based on an approach 
developed previously by Kerman et aZ.E3 Section [II summarizes the relationship of this GRPA with diagrammatic 
perturbation theory. In section IV we apply our formalism to a bilayer system at flUing factor v = 1. We evaluate the 
effect of fluctuations on the quasiparticle spectral functions and renormalization of the order parameter due to these 
fluctuations. Section ^ summarizes an adiabatic approximation for the quasiparticle self-energy in a quantum Hall 
ferromagnet and its close relationship with the two more general approaches we discuss. We conclude the paper with 
a brief discussion in section 



II. AUXILIARY FIELD FUNCTIONAL INTEGRAL APPROACH 



In bilayer systems, the anisotropy of the Coulomb interaction in pseudospin space makes it necessary to treat 
fluctuations in the Hartree and exchange channels on an equal footing. A general way to treat fluctuations around an 
ordered mean-field is the Hubbard Stratonovich transformationBil. In this approach, by introducing an integral over 
a Bose field 0, we convert the two-body fermionic interaction term into a one-body term coupled to this field. Since 
the effective Hamiltonian obtained is quadratic in fermion operators, the trace over these degrees of freedom can be 
performed exactly and we obtain an effective action for the auxiliary Bose field. Stationary phase approximations to the 
bosonic action give various mean-field theories. The nature of the mean-field state depends upon the way the Hubbard 
Stratonovich transformation is implemented and the predictions for collective excitations depend qualitatively on the 
mean-field around which Gaussian fluctuations in the Bose field are considered. For example, fluctuations in the 
exchange mean-field produce spin-wave excitations in ferromagnets, whereas fluctuations in the Hartree mean-field 
lead to plasmons in an electron gas. For a bilayer system, we need a method which treats both Hartree and exchange 
fluctuations on an equal footing. It is known that the standard HS transformation does not yield a Hartree-Fock 
mean-field and cannot capture both Hartree and the exchange fluctuationsO. 

Kerman—Levit and Troudet have presented a generalization of the HS transformation which overcomes these 
limitationsEj. To establish our notatiorO (which differs from that of Kerman et al.) we briefly review the previous 
work which describes fluctuation corrections to mean-fleld approximations for the grand potential. Our approach, on 
the other hand, enables us to systematically improve upon mean-fleld approximations for correlation functions. 

Consider a many fermion Hamiltonian in second quantized form 

H = Y^ K^pcUp + ]^Y. {o^mmcUicsc^, (2.1) 

ap aP^S 

where K is a. one-body (kinetic) term and F is a two-body interaction. Here a denotes (a set of) single-particle 
quantum numbers. For example, in the case of a bilayer system a — {n, k, a) where n is the Landau level index, k is 
the intra-Landau level index and a is the pseudospin label. Introducing pair-labels a = (a'a) and the density matrix 
Pa — cl^iCa, the kinetic term becomes KaPa where summation over repeated index is assumed. Similarly the interaction 
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term can be written as Vab^ab where (lab 



in an obvious matrix notation, we write ( 2.1) as H 
ensemble is given by 



: PaPb '■ is the normal ordered two-body term and Vab 
" Kp + Vn/2 
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(2.2) 



where T denotes the imaginary-time ordered product, N is the number operator, e = P/M and in the thermodynamic 
limit a/ (3 ^ p,, the chemical potential. Since the number operator commutes with the Hamiltonian, we can expand 

only the e^^^ term. The central idea of this generalized Hubbard transformation is to expand around an arbitrary 
two-body interaction U, thus treating Kp+ UCt/2 as the dominant term and {V — U)Cl/2 as a perturbation. Using 
the Gaussian identity 
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(2.3) 



at each time index j, a typical term in Eq.( 2.2) becomes 
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(2.4) 



Here Af^ = Af x J\f stands for a repeat sum over the Af single-particle state labels in the Hilbert space and Vcfij 
implies both a product over single-particle labels and the relevant normalization factors. Introducing these fields at 
each time-step j we get 
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(2.5) 
(2.6) 



We stress that since {(jy^) « 1/e, the limit e ^ in Eq.( 2^) must be taken after all integrals are done. 

Systematic approximations to the grand potential are obtained by expanding the bosonic action ( |2.6| ) around its 
minimum. In the limit as e — s- 0, the configuration of fields (jP which minimizes the action {dS/d(j)\^o = 0) is given by 



,Af; 



(2.7) 



where O^o indicates thermal average with the mean-field Hamiltonian hi = {K — (j>i)p. We emphasize that Eq.( 2.4) 
is an exact identity; since the variance of the auxiliary field is proportional to 1/e, the interaction term with prefactor 
is of the same order as the kinetic term with the prefactor e after all the (/)-integrals are done. The stationary phase 
approximation, however, corresponds to replacing the measure I?</)j by T)(j)jS{(j)j — (j)^). Therefore, in the stationary 
phase approximation, interaction term with the prefactor can be neglected and does not contribute to the mean- 
field solution. It is clear from Eqs.( 2.7) that the mean-field Hamiltonian is solely determined by the trial interaction 
U and is independent of the actual two-body interaction V. We get the Hartree (exchange) mean-field by choosing 
U ~ —V {U = +V'^'^), whereas —U = {V — V^) = (the antisymmetrized interaction) gives the Hartree-Fock 
mean-field. 

We improve upon the stationary phase approximation by considering quadratic fluctuations around the me an fi eld 
. These fluctuations, of course, depend upon the true microscopic interaction V . Expanding the action ( 2.6) to 
second order in the fluctuating flelds 0j = 0" + (50j gives 



eM{j,j') = 6 S,,,U-' - e (1 - %0 + ^Sj^S 



d(j)jd(l)j' 

where only time indices are explicitly shown. The matrices D and S are defined by 
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and 



{Paf3ij))4>o{pjS{j'))4,o, 



(2.10) 



where a,P are the single-particle labels and stand for the time indices. The particle-hole response function D is 
the same as the Hartree-Fock mean- field susceptibility and the fluctuation matr ix is essentially the propagator 
for collective modes around the mean-field. It is clear from Eqs.( ^Tsl ), ( 2.£) and ( 2.10| ) that properties of the collective 
modes are determined by the mean-field (jfi and the true microscopic interaction V. We note that among all trial 
potentials, only the antisy mme trized interaction optimizes the quasiparticle grand potential, flgp{U) = S[(p^] +tTS/2, 
obtained from the action ( |]^). Therefore, we will concentrate on the Hartree-Fock mean field which corresponds to 
U = -V^. 

We now extend this approach to the evaluation of the one-particle Green's function Gab{Ta,0) = +{Tca{Ta)cl{0)) . 
Starting with this definition and introducing auxiliary fields at each time step, we get 



GabiTa) = lim — Tr r 



M . s 



lim - /P0e 



(2.11) 
(2.12) 



where S is the bosonic action ( |2.6| ) and {{ca{Ta)cl{0)))^ represents the finite-e expression for thermal average with 

the auxiliary Hamiltonian hj = (K — 4)^)p. 

Systematic approximations to the exact Green's function are obtained by expanding the action S and the (p- 
dependent thermal average around the mean-field value (j)^ . At the mean-field level, we obtain the Green's function 
of non-interacting quasiparticles in the self-consistent field (jP 



(2.13) 



where w„ = 7r( 2n -|- l)//3 is a fermionic Matsubara frequency. To consider the effect of fiuctuations, we expand 



the action ( |2.6| ) to second order in fluctuating flelds 
(/)-dependent Green's function 



+ 



replace the mean-field Green's function by a 
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and expand (2.14) in powers of dcj) to all orders. The quadratic expansion of the action ( 2.6) gives non-interacting 
collective modes around the mean-field (jP. Expanding the ^-dependent Green's function gives an effective interaction 
between the mean-field quasiparticles and the collective modes. The resulting integrals over the fluctuating flelds are 
performed using Wick's theorem. In this approximation the Green's function has a self-energy that is the sum of all 1- 
particle irreducible diagrams obtained from mean-fleld Green's function C/" and the collective-mode propagator M~^. 
We approximate this self-energy by the flrst irreducible diagram, which describes a single scattering of a collective 
mode and a mean-fleld quasiparticle (Figure]^). The contribution from such a diagram is given by 



^abiiUJn) = 51 51 - «^«)M^J, j,^,(ir2„) 



(2.15) 
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where fin — 27rn//3 is a bosonic Matsubara frequency. Eq.( 2.15 ) is a general expression for self-energy of a quasiparticle 
because of its interactions with the collective modes around any mean-field. We emphasize once again that this 
calculation is to be done with a finite value of e and that the limit e — *■ is taken only after all the 0-integrals are 
evaluated. Since the dimension of the matrix M is of order 1/e and since the matrix elements of M are dependent 
on e, we have to be careful about elements of the inverse matrix. As we are interested in fluctuations around the 
Hartree-Fock mean fleld, we use U = —V^ and Eq.( 2.8) for the collective-mode propagator, and arrive at the following 
expression for the quasiparticle self-energy 



Yiabii^n 
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(2.16) 



Here we have used the identity A^j'^ = dlndct A/dAj^, and the fact that IndetAf = lndet[l -I- V^D] - tT{V^D). 
We will call ( ^.16 ) as the generalized random phase approximation (GRPA) for quasiparticle self-energy. It is clear 
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from Eq.( that the GRPA self-energy includes diagrams with n > 2 interaction terms, and in particular, it 

does not contain the mean-field self-energy contribution. This approximation for the self-energy arises naturally in 
our auxiliary field functional integral approach. In the next section we discuss the corresponding approximation in 
diagrammatic perturbation theory. 



III. GENERALIZED RANDOM PHASE APPROXIMATION 

In this section we describe the generalized ran dom phase approximation for the quasiparticlc self-energy in dia- 
grammatics language and relate it to exprfission ( 2.16| ) derived in the preceding section. In the case of particle-hole 
response functionsO or the grand potentialEj, appropriate conserving generalizations of random phase approximation 
which take into account fluctuations in both, direct and exchange, particle-hole channels have been discussed in the 
literature. 

Similar generalizations for the quasiparticlc self-energy are less than transparent. Several approaches and approx- 
imations have been used in the past. Some start from the "bubble" RPA which describes the screening and use 
vertex-corrected bubbles, while others start from the "ladders" RPA which captures the spin-wave dynamics and 
include the Hartree corrections by using screened interaction. We are not aware of a systematic generalization for the 
self-energy. 

The appropriate generalization of the RPA self-energy which we propose is summarized diagrammatically in Fig. |^. 
The physical content of this self-energy is determined solely by the nature of fluctuations encoded in the four-point 
vertex r^^^(g', We use the antisymmetrized interaction as the bare four-point vertex, Pq'*'' = V^. The self-energy 

obtained from the bare four-point vertex by contracting the incoming and outgoing labels on the top is the same as 
the the Hartree-Fock mean-field approximation for the self-energy. If instead of the antisymmetrized interaction, we 
use the direct (exchange) interaction as the bare vertex, we will get the Hartree (exchange) mean-field self-energy. We 
encode the fluctuations in direct and exchange channels into the four-point vertex by summing all the particle-hole 
ladders with the antisymmetrized interaction. In a matrix notation, the equation for the four-point vertex is 



r(4) = v'^ -v^ ■ D ■ p^"*) 



(3.1) 



where D is t he H artree-Fock mean-field susceptibility and we have suppressed all single-particle and imaginary-time 
indices. Eq.( ^j|) can be formally solved 



r(4) 



1 + V^-^ • D) 



(3.2) 



In general, the matrix (1 -f V^D) is not diagonal its indices and therefore it is not possible to obtain an analytic 
expression for the four-point vertex. The GRPA self-energy is obtained from the four-point vertex by contracting its 
incoming and outgoing lines on the top and subtracting the mean-field self-energy contribution from the resulting 
self-energy. The formal expression for fluctuation self-energy is given by 
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aa',bf3' 



(3.3) 



This is the diagrammatic equivalent of the self-energy approximation ( ^.15 ) discussed in the previous section. 

The diagrammatic content of this approximation is shown in Fig. H The ficat set of diagrams on the left, (a), 
has been traditionally used to include the effect of screening in RPA self-energyli3. This set of diagrams is generated 
by using the direct interaction as the bare four-point vertex in Eq.( ^3.l| ). The diagrams in this set can be formally 
summed and their contribution to the quasiparticlc self-energy is given by 



^ in„ a'f3' 



V 



(3.4) 



saed Coulomb interaction. The second set, (b), has been used to 
3. This set is generated by using the exchange interaction as the 



where we identify = (1 + VD) V as the sc 
include the effect of spin-waves in a ferromagnetl 
bare vertex. As expected, it represents only exchange fluctuations, which determine the spin-wave dynamics. The 
diagrams in this set, when formally summed, give 
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(3.5) 
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Eq. ( |3.5| ) reproduces the results for electron self-energy because of spin- waves in a single- layer quantum Hall system at 
ly — 1, where the physics of spin- waves is determined solely by exchange fluctuationsEj. The generalized random phase 
approximation introduced here includes another class of diagrams, (c), which represents combinations of (competing) 
direct and exchange fluctuations. We note that the set of diagrams in Fig. ^ arises naturally in the functional integral 
approach as do the subsets (a) and (b), when considering fluctuations around the Hartree or the exchange mean- field 
respectively. We refer to the sum of diagrams in all three sets as the GRPA self-energy. This approximation for the 
self-energy can also be derived by differentiating the grand potential with respect to the meaitfield Green's function 
when the appropriate generalized random phase approximation for the grand potential is usedEj. 



IV. APPLICATION TO DOUBLE LAYER SYSTEMS 



In this section, we apply the formalism developed in preceding sections to a double-layer quantum Hall system 
at total filling factor v — 1. We consider a completely spin-polarized ground state so that the layer index and the 
intra-Landau level index are the only dynamical degrees of freedom. The properties a£ihe Landau level wavefunctions 
allow us to associate a two-dimensional momentum with a pair of orbit-center labelaHil. The existence of this unitary 
transformation between the 2D wavevectors and pairs of orbit-center labels permits the significant progress that can 
be achieved analytically in the following calculations. 

The Hamiltonian for a bilayer system with interlayer tunneling amplitude At is 

Ho-flN = Y,^AaCka, (4.1) 

V^Vo + V,, (4.2) 



2 

Here, a =t= — o' denotes a symmetric state (S) while a =| denotes the antisymmetric state (AS), and e^- = —a At/2 are 
the bare single-particle energies measured from the chemical potential. The ki are the angular momenta which are good 
quantum numbers of the single-particle Hamiltonian in the symmetric gauge A = {—By/2, Bx/2, 0); Vq — {Va + Ve)/2 
and Vx = {Va — Ve)/2 are sums and differences of interlayer and intralayer Coulomb interactions. In the present 
choice of quantization axis in pseudospin-space, the interaction Vx reverses pseudospins of the scattering particles and 
therefore it does not commute with the pseudospin polarization operator. 

We start with results for the Hartree-Fock mean-field approximation. The Matsubara Green's function for the 
one-body Hamiltonian Hq is 

C/fco.(iw„) = — . (4.3) 

ILUn - Ecr 

To obtain the Hartree-Fock mean-field Green's function, we solve the Dyson equation with Hartree-Fock self-energy 

^ = - = - [ro(0)n;^(e^) + r,(0)n^^(e^)] , (4.4) 

where (A = 0, x, A, E) 

= ^ E VA(p)e-^'''/^e'^-(^~x^^)'' (4.5) 
p 

and npix) is the Fermi occupation number. The Y\{q) are interactions between an electron and an exchange-holc0 
whose center is separated by a distance qf'. The Hartree-Fock self-energy is diagonal in the orbit center and pseudospin 
labels, and the Hartree-Fock mean-field Green's function is given by 

5L(*^n) = (4.6) 



At low temperatures the stable solution of Eq.( 4.4) is given by = — [At -f m2r£;(0)] /2 = —Asas/2. Particle-hole 
symmetry ai v = 1 implies that = —£5. The exchange enhancement of the symmetric-antisymmetric energy 
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splitting, Asb = mzTE{0), is determined by the interlayer interaction and survives in the Umit of vanishing interlayer 
tunneling, giving rise to spontaneous phase coherence. The spectral function for the Hartree-Fock Green's function is 
given by Aa{uj) = 27r5(w — e°). In a mean-field approximation, the ground state is fully pseudospin polarized along the 
positive z-direction. The low-lying collective modes around this mean-field involve slow variations in the pseudospin 
field i. e. the pseudospin waves. Fluctuations out of the y — z plane in pseudospin-space correspond to the transfer 
of charge from one layer to another whereas pseudospin fluctuations in the y ~ z plane correspond to variations in the 
relative phase between states localized in the top and the bottom layer. 

We now evaluate the contribution of these pseudospin waves to the quasiparticle self-energy and its effect on the 
mean- field spectral function. The calculation of four point vertex F*^^^ or the fluctuation matrix M is particularly 
simple in the present case. Since momentum and frequency are good quantum numbers, the fluctuation matrix is 
diagonal in these indices and effectively has only the pseudospin labels. In other words, the matrix equation for the 
four-point vertex having orbital and pseudospin labels, Eq.( |3.l| ), reduces to an algebraic equation with matrices having 
only pseudospin labels. Furthermore, because of the constraints which ensure that the Hartree-Fock susceptibility D 
is nonzero, only a 2 x 2 submatrix of the entire 4x4 matrix V"^D is nonzero. This non-vanishing 2x2 submatrix is 
given by 



{aia2\V^{q}D{in„)\ 



(«.-r.)- 1 



-{vx ~ Fo' 



(4.7) 



Here 2TTpVx{q) = e"'' ' ^^Vx{q) represents the effect of electrostatic fluctuations whereas F(g) represent the effect of 
exchange fluctuations. The entire 4x4 fluctuation matrix is a trivial extension of Eq.( [4.7[ ). The roots of the equation 
det [l -|- V"^{q)D{in)\ = give the bilayer pseudospin- wave dispersionflfl iQ = Esw{q) where 



EL{q) - a{q) ■ b{q)= [AsAS - ^Aiq) + 2t;,(g)] • [AsAS - rem 



(4.8) 



Thus, the functional integral approach reproduces results obtained earlier by diagrammaticsQ and by the single-mode 
approximationlll. Figure ^ shows typical pseudospin-wave dispersions. For a finite interlayer tunneling. At ^ 0, the 
collective mode is gapped at zero wavevector because the U(l) symmetry in the y — z plane in pseudospin-space is 
explicitly broken. As d increases, the minimum in the pseudospin-wave spectrum near « 1 reaches zero at a critical 
layer separation dcr- In Hartree-Fock theory, at this point, the ground state changes from a uniformly pseudospin- 
polarized state to a pseudospin-density W9f^ statecS. This softening of pseudospin-wave at a finite wavevector has 
been associated with the phase transitionE tllj f rom a phase-coherent quantum Hall state to a compressible state. 

Using the mean-field Green's function ( 4^) and the 4x4 fiuctuation matrix we arrive at the following analytical 
expression for zero-temperature symmetric state self-energy 



Ss(«w„) 



2TrP (Es^jp) + AsA s) 



•)2 r 



e{p} - Esniip) 



Esw{p) 



^AS 



'^As{-i^n), 



(4.9) 



where e(p) — [a(p) + b{p)] /2. This remarkably simple expression is the first principal result of this work. The Dyson 
equation relating the full Green's function to the mean-field Green's function and the fiuctuation self-energy is 



(4.10) 



We obtain the quasiparticle spectral function from the retarded GRPA self-energy by analytically continuing ( 4.£) 
to real frequencies, iujn ^ to + irj, 



Asiuj) = -2 Im Gsiuj + irj) 



'2 Im 



1 



AAs{-t^)- 



(4.11) 



This spectral function has a (5-function contribution at frequency lo* which satisfies the Dyson equation cj* — = 
Y,s{lo* -f iif). The continuum piece of As{uj) is nonzero in the region where the retarded self-energy Ss(w -I- iTy) has a 
branch cut, i. e., at frequencies in the interval Is = [aj<s, u;>s] = [e^is + -E'™", e^is + As^is] . Here i?™" denotes the 
minimum of the pseudospin-wave energy near ql k 1 and As^s = Eswil — > oo) is the mean- field splitting. Explicitly, 
we can write the spectral function as 



As{uj) = 27rz5(5(w — to*) + 9{uj — (jj<)6'((jJ> — lo) 



-2 ImY^siio) 



[oo~4~ ReS]s(t^)] + [liaY.s{Lo)Y 



(4.12) 
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where = (1 — dTj s I dijj\ui-')~^ < 1 is the quasiparticle renormahzation factor. The mean-field spectral function 
As{uj) = 2t:5{uj — e^) is recovered by neglecting the fluctuation self-energy, which implies w* = and zg = 1 = zas- 
When we include the fluctuation self-energy, the spectral weight at the 5-function peak is reduced to zs < 1 and is 
distributed into a continuum piece at positive energies, cu > fi. Figure |^ shows the evolution of continuum part of 
the symmetric-state spectral function with increasing layer separation. We see that the spectral weight at positive 
energies, an indicator of the strength of quantum fluctuations in the pseudospin polarized mean-field ground state, 
increases with increasing layer separation and as required by the sum-rule, the quasiparticle renormahzation factor 
zs decreases. The presence of the spectral weight at positive energies has been detected at low temperatures by 
Manfra et al. in optical absorption experimentaHj. The redistribution of the spectral weight indicates the mixing of 
many-body states having antisymmetric quasiparticles into the completely pseudospin-polarized ground state. One 
measure of this mixing is the suppression of the order parameter 

/oo 
— npiuj) [As{uj) - AAsi^^)] (4-13) 
-oo 

from its mean-field value — 1. Using the particle-hole symmetry at = 1 we get As{lo) — Aas{~''-^)- This leads 
to the following simple expression for the zero-temperature renormalized order parameter 

m,{T = 0) = 2zs - I, (4.14) 

2nP ^ {E,^ (p) + As AS? e(p) - E,^ (p) 



zs 



1 



(4.15) 



Eq. ( |4.15| ) is the second principal result of this work. Figure ^ shows contours of renormalized pseudospin polarization 
in the parameter space (d, A^). In the mean-field approximation, the pseudospin polarization is not susceptible to the 
softening of the collective mode at a finite wavevector as d ^ dcr- When we include the effect of collective modes, 
the order parameter is strongly suppressed as the layer separation approaches the critical layer separation. We plot 
the typical dependence of renormalized polarization on interlayer spacing d in Fig. We interpret the strong 
suppression of thej-,q>£ijier parameter because of the quantum fluctuations as a precursor of the phase-transition to a 
compressible stati 



V. ADIABATIC APPROXIMATION 



In the preceding sections, we presented the generalized random phase approximation for the quasiparticle self- 
energy using a functional integral approach and discussed its equivalent approximation in diagmmmatic perturbation 
theory. In this section, we present an adiabatic approximation for the fluctuation self-energyEIl. This approximation 
is intuitively appealing and is applicable to a laxge class of problems, including metallic ferromagnets, where the 
effects of dispersive band structure are importantEZl. We will show that for bilayer systems, the fluctuation self-energy 
calculated using the adiabatic approximation is identical to the GRPA self-energy. 

The starting point for this approximation is the Hartree-Fock mean-fleld-theory Hamiltonian where we assume 
that the ground state is pseudospin polarized. The mean-field Hamiltonian couples the pseudospin polarization ruz 
to the self-consistent field (jP generated by other electrons. In the adiabatic approximation, we treat this field 4> as 
a dynamical variable and consider small fluctuations of this field from its self-consistent value cjfi. We expand the 
Hamiltonian, considered as a function of quasiparticles and the field, around the self-consistent value of the fleld. 
The leading-order term in this expansion gives an effective interaction between the mean-field quasiparticles and 
the collective modes around that mean-field. We approximate the fiuctuation self-energy by a single-coUective-mode 
exchange diagram similar to the approximation commonly used to treat phonon-exchange in metalalS. 

It is particularly easy to evaluate this self-energy for a bilayer quantum Hall system becaiise the Landau level 
bands have no dispersion and because particle-hole pair-momentum is a good quantum numberEil. We start with the 
Hartree-Fock Hamiltonian in the pseudospin polarized state 

Hhf ^ -^^crcl^Ckcr + {k'a'\UHF\ka)cl,^,Cka, (5.1) 

k,(T k'cr'ka 

where the mean-field potential UHFio) is given by 
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UHF(.q) 



^ A 




p9 T^9 



(5.2) 



Here are the Pauh matrices acting in the pseudospin-space, v\ and T\ are the direct and exchange interactions, 
and 



k'k (j' a 



"iq-r 



\k)clu 



(5.3) 



are the charge and pseudospin densities. To obtain the dyn amics of the pseudospin waves and their effective interaction 
with the mean-field quasiparticles, we linearize Eq.( 5^) around the pseudospin polarized state given by (m|.) = 



A5q^Q/2iTp . The transverse fluctuations correspond to charge- imbalance modulations with wavevector q whereas 
the fluctuations in nvt represent the variation of relative phase between states localized in the top and the bottom 



layer. These fluctuations are canonically conjugate, [TO^,m^] 
using bosonic creation and annihilation operators 



2iA/2TTp, and therefore we quantize them 



A 
2^ 



(5.4) 



The Hamiltonian for the bosonic modes introduced in ( 5.4) is obtained from the Hartree-Fock energy functional 

2nP 



Ehf [m\ = — ^"1-^0 - 



A 



■E 



yP 
^ E 



yP 

r T E 



V V 



(5.5) 



We expand energy functional ( 5.5) around the mean-field state (m^) — Adqfi/2Trl and quantize the fluctuations 



using Eq.( 5.4). A simple calculation gives the following Hamiltonian for the pseudospin-waves 

A,7 



E 



(4 



p ~p 



(5.6) 



where e{p) — [Asas + Vx — Tq] — [a{p) + b{p)]/2 and X{p) — [vx — Tx] — [a(p) — b{p)]/2. The anisotropy of Coulomb 
interaction in the pseudospin-space is reflected in the fact that the number operator = ataf does not commute with 
the Hamiltonian ( 5.6). The Hamiltonian is diagonalized by Bogoliubov transformatior@ and we get non-interacting 



pseudospin-waves with the dispersion E"^^ (q) — [if) — \^ {q) = a{q)- b{q) obtained in preceding sections. The effective 
interaction between pseudospin waves and the mean-field quasiparticles is obtained by expanding the Hartree-Fock 
Hamiltonian ( 5.1) in terms of pseudospin-wave operators. At leading order we get 



^k'k 



E^ 

p 

E^ 

p 



(4. 



sCkAS 



"i'As'^ks ■ alp-) Nk'k{p), 



(5.7) 



where the interaction matrix elements are given by 

Mk>k{q) - i{2^l^){k\e-''!-'\k') [vx{q) - T^iq)] , (5.8) 
7Vfc,fe(q) = i{2Tif){k\e-'^-^k') [vx{q) - Tx{q)] ■ (5.9) 

A straightforward albeit lengthy calculation that follows line similar to the earlier wor gives fluctuation self- 
energy expression ( identical to the one obtained using the functional integral approach or the diagrammatic 
approximation. 

VI. SUMMARY AND DISCUSSION 



This paper discusses an approximation for the self-energy of interacting electrons that goes beyond the venerable 
random phase approximation (RPA) . Depending on context the RPA describes in an approximate way the interactions 
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of electrons with fluctuations in either direct or exchange particle-hole channels. There are many circumstances, 
however, where both channels must be treated on an equal footing if the qualitative physics is to be captured. The 
objective of the approximation scheme discussed here is to fill this need. We have developed our approximation 
in the language of a path- integral approach developed by Kerman et at, that, unlike standard auxiliary-field path 
integral approaches, enables expansions to be made around Hartree-Fock mean-field states in which the quasiparticles 
experience both direct and exchange fluctuating potentials. The same approximation for the self-energy can be 
derived by differentiating the approximation for the grand potential that follows from this generalized RPA scheme, 
with respect to the mean-field Green's function. We have identified the infinite subset of many-body perturbation 
theory Feynman diagrams that are included in this self-energy approximation. As far as we are aware, this self-energy 
approximation has not been considered previously. 

We have applied our formalism to the ground state of bilayer quantum Hall systems at filling factor 1^ = 1. This 
ground state has a broken symmetry in which phase-coherence is spontaneously established between the two layers in 
order to take advantage of interlayer exchange and improve inter layer correlations. Its collective excitations involve 
both fluctuations in the interlayer phase and fluctuations in the difference between charge densities in each layer. In 
order to describe the energetics of these fluctuations, it is necessary to include both direct and exchange particle-hole 
interactions. Bilayer quantum Hall systems therefore provide an example of a situation where the standard RPA 
will fail quite badly. Our approximate self-energy properly describes the interactions between quasiparticles in this 
system and the collective excitations of its order parameter field. Bilayer quantum Hall systems provide an interesting 
example for our approach because the importance of fluctuations for the ground state properties can be tuned between 
zero and large values by varying the layer separation. It is now-weli-jestablished by both theory and experiment that 
bilayer quantum Hall systems have a quantum phase transitiorfl'M'LJ at a finite value of d where the phase-coherence 
is lost. Our approximate self-energy calculation correctly obtains an order parameter that vanishes beyond a critical 
value of d, however, we do not expect it to be reliable very close to the transition point, whether it is first order or 
continuous. The bilayer quantum Hall case is also favorable because particle-hole excitations of the ground state have 
momentum as a good quantum number and there is only one excitation at each momentum. In effect, there is only 
one excited state at each momentum, the collective excitation. Much as in the case of a one-dimensional electron 
system, there is no particle-hole continuum in the excitation spectrum. This greatly simplifies the calculation and 
allows us to derive a remarkably simple analytic expression for the self-energy. 

We have confirmed the physical content of our theory by demonstrating its equivalence, for the case of bilayer quan- 
tum Hall systems, to an adiabatic approximation in which the Hartree-Fock single-particle Hamiltonian is expressed 
in terms of the order-parameter field, via the density matrix. The possibility of specifying the density matrix that ap- 
pears in the Hartree-Fock single-particle Hamiltonian in terms of the order parameter field is unique to quantum Hall 
systems. Therefore this simple relationship cannot always be established. Nevertheless, this example demonstrates 
the ability of the approximation we discuss to accurately capture the effect of important fluctuations on quasiparticle 
properties. 

We believe that our approximate self-energy will be useful for other situations as well. A potentially interesting 
example is provided by the case of superconductors for which coupled fluctuations in the phase of the order parameter 
and of the charge density form elementary excitations and both pairing and electrostatic energies are important. In 
the case of superconductors it is particle-particle channel fluctuations are important for phase- fluctuations, unlike the 
case of bilayer quantum Hall systems. The problems can be made quite similar, however, by performing a particle- 
hole transformation for, say, the down-spins of the superconductor, turning the superconducting order into easy-plane 
ferromagnetism. In fact, going in the opposite direction, the order of a bilayer quantum Hall ferromagnet can be 
regarded as an electron-hole pair condensate by making a particle- hole transformation in one of the layers. In the 
case of a superconductor it is well-known that the long-range of the electron-electron interaction turns the collective 
modes into gapped plasmons, usually leading to high accuracy for the mean-fleld treatment of a superconductor. This 
can change, however, if the collective excitations become soft at large wavevectors as they do in bilayer quantum Hall 
systems. Such softness could be associated, for example, with partial Fermi surface nesting. 

In this work, we have concentrated on quantum fluctuation effects at zero temperature. It is easy to obtain 
corresponding expressions for finite temperature self-energy, though they are somewhat cumbersome. In the limit of 
vanishing layer separation, d 0, the finite-temperature expressions reproduce, with the appropriate identification 
of pseudospin with spin, the spin- wave contribution to the electron self-energy in a single-layer systemc£l. We have 
neglected the form-factors which encode the effects of finite well-widths since their inclusion does not change the 
results qualitatively. They can be incorporated without difficulty when modeling a specific system. 
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FIG. 1. Schematic representation of self-energy approximation ( [2.14| ). The directed solid line represents the mean-field 
Green's function Q''{iuj — iQ,) and the wavy line represents the collective-mode propagator M~^(iQ). The matrix elements at 
the vertex do no depend upon the energy transfer fi. 
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FIG. 2. Diagrammatic summary of the GRPA for self-energy. We use the antisymmetrized interaction as the bare scattering 
vertex, Fq*^ = V"^, to take into account fluctuations in the direct and the exchange channels. The self-energy is obtained from 
the scattering vertex by contracting the incoming and outgoing labels on the top. 
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FIG. 3. Diagrammatic content of the GRPA self-energy ( 2.14). The first set of diagrams (a) has been traditionally used to 
include the effect of screening wherea s the second set (b) has been used as the RPA self-energy in itinerant ferromagnets to 
include the effect of spin-waves. Eq.( 2.14) includes another class of diagrams, shown in (c), which captures the contribution 
of competing Hartree and exchange fluctuations. We refer to the sum of diagrams in all three sets as the GRPA self-energy. 
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FIG. 4. Typical collective-mode dispersion for a pseudospin polarized state. For a finite interlayer tunneling, At 7^ 0, the 
collective mode is gapped since the U(l) symmetry in the y — z plane is explicitly broken. The minimum in the pseudospin- wave 
energy near ~ 1 is because of competing Hartree and exchange fluctuations present in the term a{q). 
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FIG. 5. Continuum part of the symmetric-state spectral function: The zero of energy is at the chemical potential. As the 
interlayer spacing increases, the spectral weight at positive frequencies increases and the spectral weight zs in the (5-function 
peak at the symmetric-state quasiparticle pole decreases. This presence of spectral weight at positive energies has been detected 
by Manfra et alH 
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FIG. 6. Renormalized polarization contours. The rapid suppression of the order parameter only close to the phase-boundary 
indicates that interactions between the collective modes can be ignored for most part of the parameter space in the 
phase-coherent regime. 
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FIG. 7. Typical renormalized polarization at T = 0. The mean-field polarization is independent of the changes in the nature 
of collective mode spectrum which occur as d approaches the critical layer separation, dcr/l ~ 2.15 for At = .15(e^/eZ). The 
renormalized polarization vanishes rapidly close to the phase-boundary because of increasing quantum fluctuations. 
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